Synchronous Behavior of Two Coupled Biological Neurons 
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We report experimental studies of synchronization phenomena in a pair of biological neurons that 
interact through naturally occurring, electrical coupling. When these neurons generate irregular 
bursts of spikes, the natural coupling synchronizes slow oscillations of membrane potential, but not 
the fast spikes. By adding artificial electrical coupling we studied transitions between synchrony 
and asynchrony in both slow oscillations and fast spikes. We discuss the dynamics of bursting and 
synchronization in living neurons with distributed functional morphology. 



PACS number(s): 87.22. Jb, 05.45. +b, 87.22.-q 

The dynamics of many neural ensembles such as cen- 
tral pattern generators (CPGs) or thalamo-cortical cir- 
cuits pose questions related to cooperative behavior of 
neurons. Individual neurons may show irregular behav- 
ior [|IJ , while ensembles of different neurons can synchro- 
nize in order to process biological information Q or to 
produce regular, rhythmical activity ||. How do the ir- 
regular neurons synchronize? How do they inhibit noise 
and intrinsic fluctuations? What parameters of the en- 
semble are responsible for such synchronization and reg- 
ularization? Answers to these and similar questions may 
be found through experiments that enable one to fol- 
low qualitatively the cooperative dynamics of neurons 
as intrinsic and synaptic parameters are varied. Despite 
their interest, these problems have not received extensive 
study. Results of such an experiment for a minimal en- 
semble of two coupled, living neurons are reported in this 
communication. 

The experiment was carried out on two electrically cou- 
pled neurons (the pyloric dilators, PD) from the pyloric 
CPG of the lobster stomatogastric ganglion [||. Individu- 
ally, these neurons can generate spiking-bursting activity 
that is irregular and seemingly chaotic. This activity pat- 
tern can be altered by injecting DC current (I\ and I2) 
into the neurons, see Fig. |l|. In parallel to their natu- 
ral coupling, we added artificial coupling by a dynamic 
current clamp device Q| . Varying these control parame- 
ters (offset current and artificial coupling), we found the 
following regimes of cooperative behavior. 

Natural coupling produces state-dependent synchro- 
nization, see Fig. g. (i) When depolarized by positive DC 
current, both neurons fire a continuous pattern of syn- 
chronized spikes (Fig. ^d). (ii) With little or no applied 
current, the neurons fire spikes in irregular bursts: now 
the slow oscillations are well synchronized while spikes 
are not (Fig. @a). Changing the magnitude and sign of 



electrical coupling restructures the cooperative dynam- 
ics. (Hi) Increasing the strength of coupling produces 
complete synchronization of both irregular slow oscilla- 
tions and fast spikes (see below) . (iv) Compensating the 
natural coupling leads to the onset of independent irregu- 
lar pulsations (Fig. |2|b). (v) With net negative coupling, 
the neurons burst in antiphase, in a regularized pattern 
(Fig ; |c). 

Figure [l] summarizes the functional geometry of this 
pair J5|. Each PD cell is a motor neuron, consisting of 
a soma, a primary neurite and a neuropilar region, and 
an axon which conducts spikes to target muscles. Within 
this extended structure there is (1) frequency-dependent 
filtering of voltage signals; and (2) spatial localization of 
active membrane currents ||. First, the neurites con- 
stitute a cable. When coupled to its partner, each neu- 
ron shows an input resistance and capacitance of order 
5MQ and respectively ||. Second, the fast sodium 
and potassium channels underlying threshold-dependent 
spike generation (action potentials: amplitude w lOOmV; 
duration « 1msec) are concentrated in the membrane of 
the axon, whereas the sodium, calcium and potassium 
channels supporting slow voltage oscillations (10-30mV, 
0.3-1.0sec) are located in the neuropil ||. Cable proper- 
ties affect the passive spread of voltage signals within the 
neuron. Slow voltage oscillations experience little atten- 
uation in spreading from the neuropil to the axon, where 
they drive bursts of spikes. However, fast spike potentials 
suffer significant low-pass filtering as they spread from 
the axon to the neuropil and soma ^|. The neuropil is 
also the site of natural electrical coupling (non-rectifying 
and moderately weak: steady-state voltage attenuation 
«0.25)@. 

In the studies of synchronization in the coupled PD 
neurons both cells were active under symmetrical DC 
current injection (I\ = I 2 = I). Long records of V\(t) 
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and V2{t) were obtained, from which we show segments 
(Fig. ^) . For different values of the parameter / we see a 
regime of bursting-spiking (/ = OnA, Fig. ||a-c), and, at 
more depolarized levels, a region of pure spiking activity 
(/ = 3nA, Fig.|d). 




FIG. 1. Schematic diagram of two coupled PD neurons. 
The stomatogastric ganglion (STG) of the California spiny 
lobster, Panultrus interruptus was removed using standard 
procedures and pinned out in a dish lined with silicone elas- 
tomer and filled with normal lobster saline [jij]. The STG 
remained connected to its associated anterior ganglia, which 
provide activating inputs J3j. Separate, glass microelectrodes 
(filled with 3M-KCI; tip resistance 10-20 Mil) were inserted 
in the soma of each neuron, for intracellular voltage record- 
ing or current injection. Each microelectrode was served by 
a separate amplifier and current source. Measured voltage 
signals were digitized at 5000 samples/sec. The two PD neu- 
rons remained coupled to each other by their natural elec- 
trical synapses, but were isolated from the rest of the CPG 
by blocking chemical input synapses with picrotoxin (7.5/iM) 
and photo-inactivating other, electrically-coupled neurons Q . 
Artificial electrical coupling was provided by injecting equal 
and opposite current I a into the two neurons, such that 
la = 9a(Vj — Vi), where g a is the added synaptic conductance 
and Vi is the membrane potential at the soma of PD^ |H| . 



regimes were also seen after isolating a single PD neuron 
by photo-inactivating its partner. 
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FIG. 2. Regimes of oscillations in two coupled neurons. 
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FIG. 3. Time series for an isolated PD for different values 
of I\. PD2 is inactivated by large negative current injection. 



Our results indicated that spiking and bursting- 
spiking regimes of activity arise from the autonomous 
dynamics of individual PD neurons. In the experiment 
(Fig. [3|a, b) we recorded the membrane potential of one 
PD cell PD\ when its partner PD2 was deactivated by 
DC hyperpolarization to — 80mV, effectively suppressing 
its neural activity. The DC current I\ injected into PD\ 
was varied. At 1\ = 2 nA the activity consisted of ape- 
riodic slow oscillations surmounted by spikes (bursting- 
spiking), see Fig. |3|a. At I\ = 5nA, the neuron gener- 
ated fast spikes alone, see Fig. ^b. Thus, the voltage- 
dependent spiking and bursting properties of a single 
PD resembled those of the active pair (although the val- 
ues of I\ are shifted relative to those of a symmetri- 
cal pair, due to the shunting action of the deactivated 
neuron) (Fig. ||a,d). Similar, voltage-dependent activity 



The records presented in Fig. ga clearly indicate the 
synchrony of bursts in the naturally coupled neurons. For 
more detailed analysis of synchronization of these aperi- 
odic bursts we adopt the technique developed for the ex- 
perimental studies of chaos synchronization in electronic 
circuits fiofl . To study the synchronization of slow bursts 
we suppress the spikes in the recorded signals using low- 
pass filter with cut-off frequency 5Hz and analyze the 
"slow trajectories" given by the filtered signals VFi{t) 
and V~F2(t)- The projections on to the planes of variables 
(V F1 (t),V F2 (t)) and (V F1 (t),V F1 (t + t d )) shown in Fig. | 
characterize the level of synchrony of bursts in the neu- 
rons and the complexity of the bursts dynamics, respec- 
tively. To quantify synchronization, we calculate the dif- 
ference V F d (t) = V F x (t) — V F 2 (t) , and study the normal- 
ized standard deviation <tn — ov FD /ov F1 and normalized 
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maximal deviation A N = \V F D \ max / {V£\ ax - Vpf n ) as a 
function of g ai see Fig. |5[ 
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FIG. 4. Phase portraits of the slow components of os- 
cillations in coupled neurons plotted in the planes vari- 
ables (V>i(t),V> 2 (t))-left, and (Vpi(t), V F i(t + t d ))-right. 
t d = 0.3sec. 1=0 
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FIG. 5. on and Ajv as a function of the conductivity 
through the electrical clamp g a . 

The dynamics of slow oscillations changed as the effec- 
tive coupling conductance was altered by adding artificial 
coupling, g a . With natural coupling (g a = OnS) the slow 
oscillations stayed synchronized (Fig. Ha) despite very 
complex dynamics (Fig. ^b) (cf. Fig. ga). Additional 
dissipative coupling (g a > OnS) increased the level of syn- 
chrony between the neurons, while compensation of natu- 
ral coupling (g a < OnS) led to desynchronization (Fig. ||). 
The desynchronized, slow oscillations remained complex 
and aperiodic (Fig. ^c,d, see also Figs. 0b). Adding fur- 
ther, negative coupling conductance (g a < — 240nS: prob- 
ably overcompensating the natural synapse) caused the 



neurons to become synchronized again, but in antiphase 
(Figs. §3 , |). This regime of antiphase synchronization 
was characterized by the onset of more regular, "almost 
periodic" bursts (Fig. [|f). 

Next we describe the synchronization of the fast, spike 
oscillations. The standard criterion for identical syn- 
chronization fails here because of small fluctuations in 
spike timing. Therefore, we applied a different analysis. 
For each membrane voltage we located the times of spike 
peaks, for PDi and t+ for PT>2, and calculated the 



intervals T; = t 



(i) 



i+1 



and tv = 



tj. Tj measures 



interspike intervals in PD\, while Tj tells us about the 
difference in spike timing in PD\ versus PZ?2- If l^l did 
not grow with time and max{|Tj|} < min{Tj}, we con- 
cluded that the neurons spike synchronously. We also 
analyzed the level of synchronization by measuring the 
phase relation between spikes. In this analysis we plot- 
ted a histogram of the phase, A$j, of the ith spike of 
PD2 within the interval formed by the neighboring pair 
of spikes in PD\ (designated jth and kth, respectively), 
using the function, A$; = 180 o (if } - tf ] )/\tf ] - t { £>\. 
The results of these analyses are shown in Fig. |(| 
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FIG. 6. Analysis of synchrony in the fast spiking activity 
measured for three different regimes of neuron activity de- 
scribed in the text. Distributions of the phase lags between 
the spiking of the neurons (left). The evolution of interspike 
intervals n and T (right). 

With natural coupling, spikes were synchronized dur- 
ing tonic firing. Figure ^ and b shows such a syn- 
chronous regime (7 = 3nA and g a = 0) where all val- 
ues of A<£> are within the interval (—40°, 130°), and the 
value of Tj oscillates but remains smaller than Tj, see 
Fig. ||d. In the spiking-bursting regime, natural coupling 
did not synchronize the spikes (unlike the slow oscilla- 
tions) (Figs. ||c and d). In this condition (I = OnA and 
g a = 0: cf. Fig. 0a), the constant drift of Tj indicated 
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a difference in spike frequencies; nevertheless, the non- 
uniform distribution of A$; indicated that the neurons 
were not far from the threshold for spike synchroniza- 
tion. Indeed, the spikes of bursting neurons became syn- 
chronized when artificial (positive) coupling was added 
(Figs. |e and f: 1 = and g a = 190nS). 

Our observations indicate that the slow oscillations and 
fast spikes of these two neurons have different thresholds 
for the onset of synchronization. This can be understood 
in terms of the different sites of origin of the two types of 
voltage signal, the different mechanisms of synchroniza- 
tion, and the different conduction pathways and atten- 
uation factors involved (cf. Fig. [I] and associated text). 
The slow voltage oscillations that underlie bursting ac- 
tivity arise as a result of voltage-dependent ion channel 
activity in the membrane of neuropilar processes. The 
summed voltage signal will suffer some attenuation as 
it spreads by local current flow in the leaky cable ar- 
ray of the neuropil. However, two factors favor its ef- 
fective transmission between the neurons: the location 
of electrical coupling sites close to the site of slow wave 
generation, and the slow timecourse of the voltage sig- 
nal itself. In combination, these should allow a relatively 
strong and continuous interaction between the irregular 
slow oscillators. This mechanism resembles the synchro- 
nization seen in dissipatively coupled chaotic electrical 
circuits In contrast, fast spike signals suffer strong 
attenuation as they spread between the spike initiation 
zone at the origin of the axon and the coupling sites in the 
neuropil. These factors argue for weak current flow be- 
tween spike generators. If the spike generator of one neu- 
ron is close enough to its threshold, the transient current 
from the coupling pathway may drive it to phase-locked 
firing. In electrical circuits, this type of chaotic pulse 
synchronization is known as threshold synchronization 
[pL l|| . With natural coupling, this threshold mechanism 
can synchronize spike activity in tonic firing but not in 
the bursting regime. When the neurons generate slow 
voltage oscillations, ion channel open in neuropilar pro- 
cesses, decreasing the membrane resistance: this shunts 
the spike-evoked currents as they flow in their coupling 
pathway, causing a failure in threshold synchronization. 

As the strength of net coupling is decreased, the slow 
oscillations remain irregular with little change in wave- 
form, but make a sharp transition from synchronous to 
asynchronous behavior, see Figs @, § When the net 
coupling reaches an expected, negative conductance, the 
slow oscillations re-synchronize in antiphase and become 
regular. These bifurcations argue for a dynamical origin 
of the irregular neuronal activity. Based on these obser- 
vations we have built a two-compartment model of the 
stomatogastric neuron. The model incorporate six ac- 
tive ionic currents distributed in soma-neuropil and axon, 
and takes into account slow, intercellular Ca ++ dynam- 
ics. Two such model neurons, when electrically coupled, 
reproduce all five types of behavior found in our experi- 



ments and the transitions between the regimes are con- 
sistent with the observations reported here, see Ref . [^2) . 
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